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1. Introduction 

In physics different theories are often formulated in different languages. For example, 
quantum mechanics is formulated via states and operators in a Hilbert space endowed 
with a commutator structure. Classical mechanics for a system with n degrees of freedom 
is instead formulated in a phase space labeled by n coordinates qi and n canonical 
momenta p- 7 , which we will collect together in a unique Lp a = {q-* ,p*), with the index 
a running from 1 to 2n. This phase space is endowed with a Poisson brackets (pb) 
structure, which can be written as: {tp a , f b } p h = u ab , where the symplectic matrix u ab 



is the following 2n x 2n antisymmetric matrix: u = ^ ^ ^ J . 

In most textbooks a theory is usually developed in only one way and the alternative 
formulations of the same theory are often not mentioned at all. All this can lead 
students to the misconception that there is a very strict link between a theory and 
its mathematical formulation. Of course, there is often a particular language which is 
more suitable than others to express a physical theory, but it is important to realize 
that such a language is often not the only possible one. In this paper we want to stress 
the distinction between the physical content of a theory and the language which is 
commonly used to formulate the theory itself. In fact, we will provide a simple example 
to show explicitly how the same language (the path integrals) can be used to formulate 
different theories, like classical and quantum mechanics. 

As a consequence, we shall also show that the well-known Feynman diagrams, 
used mostly in quantum field theory to develop perturbation theory, can also be used 
in classical mechanics to provide a perturbative solution to the classical equations of 
motion. This paper can then be useful as an introduction to the main ideas and concepts 
of Feynman diagrammatics, which are studied in quantum field theory courses. 

2. A brief introduction to a path integral approach for classical mechanics 

As we said in the Introduction, the main purpose of this paper is to implement 
perturbation theory in classical mechanics via the same tools that are used in quantum 
field theory. To reach this aim we have first of all to introduce also in classical mechanics 
concepts like path integrals and generating functionals. 

It is well-known that path integrals were first introduced in quantum mechanics by 
R.P. Feynman [1]. He discovered that the transition amplitude of going from a certain 
point qi at the initial time to another point qj at the final time T can be written as 
a sum over the paths joining these two points on an interval of time T. In particular, 
it is possible to show that all the paths in phase space give a contribution to such a 
transition amplitude according to the following formula: 



where the double prime over V indicates that the functional integration is over functions 





■T 



dtL 



(1) 
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q(t) with fixed end points, i.e. g(0) = and q{T) = qj, and L is the following 
Lagrangian: 

2 

L = pq-H, H = P - + V{q). (2) 

Let us now pass to the path integral formulation of classical mechanics developed 
at the end of the 1980s The starting point of this formulation is the following: which 
is the probability density P(<pf,<pi]T) of going from a certain point of the phase space 
ifi at time to another point iff at time T7 Of course, we know that the answer is very 
different from the quantum case. In fact, in classical mechanics, once we fix the initial 
conditions (p(0) = <pi, there is only one path in phase space which solves the equations 
of motion, let us indicate it with (p^(t; (pi), where "cl" stands for "classical path". As 
a consequence, the probability density P(<ff, tpf, T) is different than zero only when the 
point iff lies on the classical path (p^(T\ (pi): 

P{<p f ,wT) = 6[<f f -<A{TM)]. (3) 

Now we want to rewrite the previous result as a path integral over all the possible 
trajectories in phase space <p(t), where the only path which contributes is the classical 
one. This can be realized via a suitable functional Dirac delta, i.e.: 

P(<p f , <p i; T) = / V"<p5[ip a - <^(t; ify)]. (4) 



Vi 



This is a first expression for the path integral for classical mechanics [2|. Let us now 
remember the following property of the Dirac deltas: if a function f(ip) has a zero 
in (p , then the following equation holds: 5(<p — <p ) = S[f((p)] ■ |/'(^ Q )|. In our case 
the function f{<p) is given by the LHS of the classical Hamilton equations of motion 
ip a — uj ab dbH = 0. As a consequence, the zero tp of f(tp) becomes the path (p^(t;(pi) 
which solves the classical equations of motion. The modulus of /' becomes instead a 
functional determinant and the weight of the classical path integral (pij can then be 
rewritten as follows: 

6& a - ¥&(t; ifi)} = 5& a ~ uJ ab d b H)det(d t 6 a b - u ac d c d b H). (5) 

It can be proved [2j that the functional determinant in the previous formula is 
independent of the phase space variables tp, so it can be omitted at this level. Using 
the Fourier transform, the Dirac delta of the equations of motion can be rewritten as a 
path integral over the auxiliary variables A a = (X q , X p ): 



5& a - u ab d b H) = j m e i fo T <x^ a -» ab d»H) 



So the path integral for classical mechanics j2] becomes a functional integral over an 
extended phase space (ip, A): 



Vf 



T 

i I dtC 



o 



(6) 



P((p f ,(pi;T)= / V'tpVXexp 

where the Lagrangian and the Hamiltonian associated with each path are given by: 

£ = X a if a - H, H = X a u ab d b H. (7) 
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In this way we have rewritten the weight of the path integral for classical mechanics 
in an exponential form, just like in the quantum case. The main differences between 
the two path integrals for quantum and classical mechanics (JTJ) and (jBJ) are in the space 
of functions which are integrated over and in the arguments of the exponential weight, 
which are defined in (J2J) and (|IJ) respectively. Nevertheless, since we have classical 
mechanics formulated in the same language of quantum mechanics, it is reasonable 
to expect that, even in classical mechanics, we can develop perturbation theory via 
Feynman diagrams. 



3. The classical generating functional 



In Eq. © we have found a path integral for the classical transition probabilities. What 
we want to do now is to define the analog of the generating functional jHj for classical 
mechanics. Let us introduce two currents, J and A, associated with q and X p respectively, 
and define: 



Z[J,A] = I d<p f H V"^V\e i ti dt{c+Jq+AX v ) . 



(8) 



When J = A = the generating functional Z[J, A] is normalized to one. In fact, using 
® and ©: 



Z[J,A) 



J,A=0 



d<p f P((p f ,<pi;T) 



The generating functional Z[J, A] can be used to formally write down the solutions of the 
equation of motion. In fact, as we will prove in Appendix A, the functional derivative of 
Z[J, A] w.r.t. J{t) reproduces exactly the solution of the equations of motion associated 
with the initial conditions y?(0) = ipf 
.SZ[J,A] 



6J(t) 



<t <T. 



(9) 



J,A=0 



Of course, for a generic potential it is impossible to calculate exactly the path integral (0) 
and consequently its derivative to get the exact expression for q c \(t; </?j). Nevertheless, 
as we will prove in Appendix B, we can derive the explicit form of the generating 
functional Zq[J, A] for a harmonic oscillator described by the Hamiltonian: 



Ha 



■ 2 2 

-u q . 



The result is: 

Zq[J, A] = exp 



dtJ{t)q {t;ipi)-i / dtdt' J(t)G R {t-t')A{t') 



(10) 



where go(^ <Pi) is the solution of the harmonic oscillator equations of motion with initial 
conditions tpi = (qi,Pi): 



Pi 

q (t; ifi) = qi cosut H — - smut 
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and is the following retarded Green function (see Appendix B for more details): 

, sin u;(t — t') 



G K (t - t') = Q(t - t' 



As a consequence, from Eq. (fTUj) we have that 
,8Z [J,A] 



6J(t) 



qo(t; <fi 



< t < r, 



(12) 



(13) 



J,A=0 



which confirms the result © in the case of a harmonic oscillator. In a similar way, we 
can get the retarded propagator (fT2"J) by acting on the generating functional with two 
functional derivatives as follows: 
, 2 S 2 Z [J,A] 



-i 



5J(t)SA(t') 



iG R (t - t') 



(14) 



J,A=0 



The generating functional (jlUJl can be considered as a starting point to implement 
perturbation theory for classical mechanics. In particular, we will call propagators 
the quantities in Eqs. (|13| ) - (jl4j) . which are coupled with the external currents in the 
generating functional. These propagators can be graphically represented according to 
the following Feynman rules: 

. = q (t;ipi) 

• . = iG R (t-t'). 



Next, let us consider a system described by a Hamiltonian H which can be split into 
the sum of the Hamiltonian Hq of the harmonic oscillator and a term of perturbation 
gV(q): 

H = H + gV(q). 
Then the TC of Eq. (JJJ) can be split as: 

H = X a u ab d b H - gX p V'(q) =H + gV(q, X p ), 

where in the first step we have used the fact that the perturbation V is independent of p. 
The generating functional associated with the system described by the new Hamiltonian 
H can then be written in terms of the free generating functional (jlUj) as follows: 



dtpf V'tpVXexp 



Pi 



i I dt(X a (p a 
o 



Hq - gV{q, X p ) + Jq + AA P 



(15) 



exp 



-w 



dtV -i 



8 



v SJ(ty 6A(t)_ 

In the last line of Eq. (fTK|) we have replaced the arguments of V, i.e. q and X p , with the 
functional derivatives w.r.t. the associated currents J and A respectively. 

According to Eq. (fTK|) . the generating functional Z can be approximately calculated 
by expanding the differential operator acting on Z up to the desired order in g. As 
we will see in the next section, this will lead to another Feynman rule which crucially 
depends on the details of the perturbation. 
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4. An example: the anharmonic oscillator 

We will now consider an anharmonic oscillator. In this case the perturbation is 

g 4 1 
V(q) = — , which implies: V(q,\ p ) = — gA p g 3 . The effect of this perturbation can 

be taken into account by introducing a new graphical element known as vertex. In 
particular, we will associate every factor in the perturbation V with an external leg: a 
dotted one, to represent X p , and three full ones, to represent each factor q. In conclusion, 
the Feynman rule for the vertex is the following one: 




Now, as a first example, let us calculate, up to the first order in g, the solution of 
the classical equation of motion for an anharmonic oscillator. First of all, let us indicate 
with a circle the generating functional Z[J, A], which has the following diagrammatic 
expression: 



Z[J,A] = 




The functional derivative w.r.t. J(t) can be indicated as: 



— i- 



. 8Z[J,A] 
SJ(t) 



J,A=0 t 




Using Eq. JSJ) we have that the solution of the classical equations of motion for an 
anharmonic oscillator can be represented via the following Feynman diagram: 



Qd(t) 



t 




+ 



t t 

which, using our Feynman rules, can be translated in the following equation: 

?ci(t) = go(*) " | f T dt' G R (t - t')q 3 (t') + 0{g 2 ). 
6 Jo 

Let us solve explicitly the previous equation, considering for simplicity the particular 
initial condition pi = 0. Then, using (|llj ) -([12 p . we have the following solution for the 
classical equation of motion of the anharmonic oscillator: 

,3 r l 



9Qj 
48u; 2 



3ut sin ut + - (cos ut — cos 3ut) +0(g 2 ). (16) 



This is the same result that can be obtained using the standard method of secular 
perturbation theory, see for example jlj. 
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5. Equivalence with the secular perturbation theory 

The equivalence of the result (fTo]) with the one that can be obtained via the secular 
perturbation theory is not an accident. In fact, it is quite easy to show that our Feynman 
diagrammatics is just a stenographic expression for the secular approach to perturbation 
theory in classical mechanics. In fact, let us consider the classical equation of motion 
for the anharmonic oscillator: 

q(t)+u 2 q{t) = - 9 -q\t). (17) 

The secular approach to perturbation theory in classical mechanics jlj consists, first of 
all, in expanding q(t) in series of g: 

q{t) = q (t) + gqi{t) + g 2 q 2 {t) + ■ ■ ■ . 

Replacing the previous expansion in the equation of motion (fTTj) we get an infinite set 
of coupled equations for the components qk(t): 

q (t)+u 2 q (t)=0, 

q\(t)+u 2 qi (t)= - l -ql(t), (18) 

q2(t)+LJ 2 q 2 (t)= - ^ql(t) qi {t). 

Then the boundary conditions g(0) = qi and q(0) = Pi can be turned into the following 
set of conditions: 

9o(0)=ft, qo(0)=Pi, 

(19) 

9fc(0) = q k {0) = Vfc>l, 
which can be used to solve the equations of motion (jl8j) . The first equation is the 
unperturbed equation, which can be solved to find the go of Eq. (JTTJ). This expression 
can be inserted into the second equation of (JTHJ). Imposing the boundary conditions 
(fT9"|. one then finds the solution for gq\{t): 

gqi(t) = - 9 -j\t'G K {t-t>)ql{t'), (20) 
which is just the contribution given by the following Feynman diagram: 




t 

This procedure can be iterated. For example, using ()2U|) in the third of Eqs. (fTHj) . we 
get: 

,2 rT 



[ T dt' G R (t - t')q 2 (t') ■ f dt" G K (t> - t")ql(t 
Jo Jo 



1 rT 

g 2 ' 



12 

which corresponds to the following Feynman diagram multiplied by 3: 
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t 




The factor 3 (called symmetry factor) is due to the fact that the vertex in t" can be 
attached to each one of the 3 continuous legs coming out from the vertex in t'. It is 
interesting to note that in this diagrammatics the only contribution to the solutions 
of the equations of motion is given by tree diagrams. This is due to the fact that the 
propagator Gr(£ — t') is proportional to smuj{t — t'), so it cannot generate loops because 
it is identically zero for t = t'. These tree diagrams are similar to the Wyld diagrams 
introduced in the theory of turbulence, see for example j3]. 

A more compact proof of the equivalence between our perturbation theory via 
Feynman diagrams and the more standard secular perturbation theory goes as follows. 
It is easy to realize that the following iterative relation holds: 



By deriving the previous equation twice w.r.t. t and using (jB.3|) . we get exactly the 
equation of motion (jl7|) of the anharmonic oscillator. This completes the proof of 
the equivalence of our diagrammatics with the standard secular approach to classical 
perturbation theory. 

6. Conclusions 

In this paper we have shown how Feynman diagrams can be used also in classical 
mechanics as a useful stenographic tool to implement the standard secular perturbation 
theory. Of course, the formalism that we have presented here can be easily modified 
to take into account other examples. In particular, a change in the free generating 
functional implies a change in the form of the propagator, while a change in the 
perturbation potential V(q) implies a change in the form of the vertex. We think 
that this approach to perturbation theory can be used to convey to students the idea 
that path integrals and Feynman diagrams are powerful tools, which can be fruitfully 
employed not only in quantum field theories but also in other areas. 




The equation associated with the previous diagram is: 
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Appendix A. 

In this Appendix we want to prove explicitly Eq. 0. Let us start from the definition 
of the generating functional for classical mechanics (JHJ): 



Z[J, A] = J dip f V"(pVX e l dt ( £ +^+ aa p) . 



Calculating the functional derivative of Z w.r.t. J{t'), and putting afterwards the 
currents J and A equal to zero, we get: 
.SZ[J,A] 



—i- 



SJ(t') 



J,A: 



= I d(pfV"(pVXq(t') Atc . 
=o J 



As usual, the functional integral over A can be performed explicitly to get the functional 
Dirac delta of the equations of motion: 
.SZ[J,A] 



SJ(t f ) 



= ! d<p f V"<pq(t') 5(if a - uo ah d b H). 

J,A=0 J 



Next, we can use the analog of Eq. (jSj) and pass from the delta of the equations of 
motion to the one of the solutions of the equations 5[(p a (t) — ip%(t] ipi)]. Also the path 
integral over (p can be performed explicitly and we get: 
.SZ[J,A] 



= q c i{t';ip i ) I d<p f 6[<p a f -&(T]<Pi)] 

J,A=0 



5J(t' 

= q cl (t'; ifi). 

This completes the proof of Eq. Q. 
Appendix B. 

In this Appendix we want to prove Eq. (fTU|) . i.e. we want to derive the generating 
functional Zq[J, A] for the harmonic oscillator. First of all, let us introduce a new 
Hamiltonian H' Q = H + Aq. With this definition the generating functional of the 
harmonic oscillator can be written as follows: 

J / T dtiXa^-Xa^dtH^+Jq) 



Z [J,A] = J diffV'ipVXj 

dip f V"ip5(ip a - uj ab d b H' Q ) e i £ AtJ \ 



where in the last step we have performed the path integral over A. Using the analog of 
Eq. (jSJ) for the Hamiltonian H' Q and disregarding the functional determinant, which is 
independent of the variables q and X p , we can rewrite the generating functional as: 



Z [J,A] =N J dip f V"y5{p a -y'°i{t- ^)}e^o 



M exp ( % dtJq' cl ) . (B.l) 
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The condition Z[J, A = 0] = 1 fixes M = 1. Finally, we have to find the explicit form of 
q' cl by solving the equations of motion associated with H' , i.e.: 

q = p, p = —u 2 q — A, 

which are equivalent to the following second order equation for q: 

q + uj 2 q = -A. (B.2) 

Let us consider the so called retarded Green function Gr(£ — t') of Eq. (fl"2*J) . which 
satisfies the following equation: 

(d 2 + uj 2 )G R (t - t') = 5(t - t') (B.3) 

with the boundary condition Gr(£ — t')=0 for t < t! . Then Eq. (1B.2|) can be solved in 
terms of the Green function as follows: 

q' cl (t) = Qo(t) - [ T dt'G R (t-t')A(t'), (B.4) 

where q (t) is the solution (fTTj) of the homogenous equation q + uj 2 q = with initial 
conditions q(0) = qi and g(0) = pi. Finally, replacing (JB.4J) into (jB.ljl . we get the 
explicit form (fTIH) of the generating functional for a harmonic oscillator. 
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